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Abstract 

We have studied the evolution of an excited electron bubble in superfluid ^He for several tens of 
picoseconds combining the dynamics of the liquid with an adiabatic evolution for the electron. The 
path followed by the excited bubble in its decay to the ground state is shown to strongly depend 
on pressure. While for pressures below 1 bar the IP excited electron bubble has allowance for 
radiatively decay to the deformed ground state, evolving then non-radiatively towards the ground 
state of the spherical electron bubble, we have found that above 1 bar two distinct baby bubbles 
appear in the course of the dynamical evolution, pointing to a different relaxation path in which 
the electron may be localized in one of the baby bubbles while the other collapses, allowing for 
a pure radiationless de-excitation. Our calculations are in agreement with experiments indicating 
that relaxed IP bubbles are only observed for pressures smaller than a critical one, of the order 
of 1 bar, and that above this value the decay of the excited bubble has to proceed differently. A 
similar analysis carried out for the 2P bubble shows that the adiabatic approximation fails at an 
early stage of its dynamical evolution due to the crossing of the 2P and IF states. 

PACS numbers: 47.55.D-,67.25.du,33.20.Kf,71.15.Mb 
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I. INTRODUCTION 



Electron bubbles (ebubbles) in liquid helium are fascinating objects with an apparently 
simple structure that have been the subject of a large number of experimental and theoretical 



studies, see e.g. Refs. 



12 



and references therein. 
The imaging of individual ebubbles moving in the liquid, some unexplained events in 
cavitation experiments,^^ and the efforts in creating and detecting multielectron bubblesi^— 
are recent issues calling for a dynamical description of the electron bubble, but not the only 
ones. For instance, the equilibration of the electron bubble in superfluid liquid helium, stud- 
ied in detail by Eloranta and Apkarian^ within time-dependent density functional theory, 
also needed of an accurate dynamical description. The ebubbles addressed in that work 
are spherically symmetric, which made the calculations affordable for using the best avail- 
able density functional (DF) for ^He, the so-called Orsay-Trento (OT) functional.^'' Other 
dynamical studies have resorted to much simpler approaches inspired on local functionals 
of the kind proposed by Stringari and Treiner long time ago,^^ or on generalizations of the 
Gross-Pitaevskii equation to the description of liquid helium.—"— They have allowed to carry 
out dynamical studies involving non-spherical ebubbles, and their interaction with vortices 
in the superfluid. However, these local approaches do not describe the superfluid accurately. 
In particular, its elemental excitations are poorly reproduced. To circumvent this short- 
coming, non-local extensions have been proposed^^ and applied e.g., to vortex nucleation in 
superfluid helium.— 

Another problem requiring a dynamical treatment, still not addressed in full detail, is 
the relaxation of an ebbuble after being excited by photoabsorption, which constitutes the 
subject matter of this work. This process couples the fairly slow displacement of the helium 
bubble with the rapid motion of the electron it hosts, producing excitations in the liquid 
that take away a sizeable part of the energy deposited in the ebubble during the absorption. 
The emission spectrum of the electron bubble after it has relaxed around the excited electron 
state has been calculated.— However, whether and how these full relaxed configurations 
are attained before decaying by photoemission was not elucidated. 

In this work we attempt a theoretical description of the evolution of the excited ebub- 
ble based on the zero temperature DF approach, using an as much accurate as technically 
feasible description of the liquid, and an electron-helium interaction that have been proved 
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to reproduce the experimental absorption energies of the ebubble. The initial configuration 
is determined by a static calculation of the excited ebubble. This state has a large radia- 
tive lifetime, of the order of several tens of microseconds, in contrast with the short time 
scale for the helium displacement, of the order of picoseconds. The subsequent dynamical 
evolution of the ebubble is described within the adiabatic approximation, which is valid for 
a period of time difficult to ascertain,^'' that we shall discuss in some detail. We will show 
that, depending on the initial excited state and the external pressure applied to the liquid, 
the bubble may keep its initial simply connected topology, or evolve towards a non-simply 
connected one made of two baby bubbles that share the probability of finding the electron 
in,-^ the electron eventually localizing in one of them while the other collapses. To reduce 
the numerical effort to a reasonable amount, we shall mostly discuss results for the collapse 
of an ebubble starting from the spherical IP state. Results for the collapse of the 2P state 
will be also shown. 

This work is organized as follows. In Sec. II we describe our model and present a 
quasi-static study of the ebubble, completing the results we have presented elsewhere, and 
recalling some technical details about the method we have used to solve the variational 
equations for the ffuid and the electron. In Sec. Ill we present the adiabatic evolution of 
the ebubble for two selected values of the liquid pressure. The validity of the adiabatic 
approximation is analyzed in Sec. IV, and a summary is presented in Sec. V. 

II. QUASI-STATIC DESCRIPTION 

We first address some properties of excited electron bubbles in liquid ^He using the Orsay- 
Trento density functional including the terms that mimic backfiow effects and are crucial to 
quantitatively reproduce the experimental phonon-roton dispersion relation in bulk liquid 
^He at zero pressure.^'' They have no inffuence on the statics of the system, and have been 
often neglected^^'^^ in the dynamics because their inclusion makes the dynamical calculations 
very cumbersome. -j^"— In practice, we have found that these terms have little effect on the 
dynamics of the electron bubble presented in this work. 

The electron-helium (e-He) interaction has been modeled by the Hartree-type local effec- 
tive potential derived by Cheng et al.— This allows us to write the energy of the electron- 
helium system as a functional of the electron wavefunction $(r) and the ^He effective 
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macroscopic wavefunction ^^(r) = y p(r) exp[z5'(r)], where p(r) is the particle density and 
v(r) = hV S{y) /rriHe is the velocity field of the superfiuid: 



2 rriHe J J 

/ dr|V$(r)p+ / dr|$|Ve_He(p) - (1) 

Irrip J J 
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In this expression, E{p) is the He 'potential' energy density, and the e-He interaction 



34. 



. In 



Ve-He{p) is written as a function of the helium density. Details are given in Refs. Iq. 
the absence of vortex lines, S is zero and E becomes a functional of p and $. Otherwise, 
one has to use the complex wavefunction ^^(r) to describe the superfiuid. 

For a given pressure (P), we have solved the Euler-Lagrange equations which result from 
the variation with respect to and $* of the zero temperature constrained grandpotential 
density $) = cl;(^^, $) - with 

= -^|V^(r)|^ + ^(p) 

2mHe 

+ ^|V$r+|$|Ve_^e(p)-MP, (2) 

z me 

where p is chemical potential of the liquid. The variation of the above functional yields two 
coupled equations that have to be selfconsistently solved 

Av^ + H + |$r^^%i^^U = ^v, (3) 



2 ruHe i h dp 



^ /\^ + Ve-He{p)^ = e^, (4) 
ZTTZe 

where s is the eigenvalue of the Schrodinger equation obeyed by the electron. 

Our method of solving the variational equations is based on a high order discretization in 
Cartesian coordinates of the differential operators entering them (13-point formulas in the 
present case) , and the use of fast Fourier transformation techniques^ to efiiciently compute 
the convolution integrals in uj{p)^ such as the mean field helium potential and the coarse- 
grained density entering the definition of the correlation energy. This allows us to use a 
large spatial mesh-step, of about 1 A size, without an apparent loss of numerical accuracy 
when we compare our results with others (see below) obtained using 3-point formulas for 
the derivatives that, as a consequence, require a rather small mesh-step to be accurate. The 



density at the boundary of the three-dimensional (3D) 140 A x 140 A x 140 A box used to 
carry out the calculations is fixed to the value of the bulk liquid density at the given P. We 
recall that knowledge of S{p) allows to determine the equation of state of the bulk liquid 
and its chemical potential, since /x = 88 /dp and P = —£{p) + pp. Eqs. (IHHl]) have been 
solved employing an imaginary time method, and we have carried out the appropriate 
tests to check the stability of the solutions. We mention that the energies we have obtained 
for the IS^IP and 1S^2P transitions^ are in very good agreement with experiment,^ and 
that our results compare well with those obtained by Eloranta and Apkarian^ using the 
same functional but a different numerical method and e-He interaction. This constitutes an 
excellent test not only for the numerics, but also for the physical ingredients employed in 
both calculations. We have recently discussed the eflFect of the presence of vortices on the 
absorption spectrum of ebubbles attached to them.— 

Upon excitation to the IP state by light absorption, the ebubble experiences a drastic 



change of shape. This is due to the fairly large radiative lifetime o 



to be 44 /xs in Ref. 26, 60 /is in Ref. Il0|, and 56 /is in Ref. l25|) as compared to any 



this state (calculated 



characteristic helium timescale, allowing the liquid to relax around the excited state. As 
a consequence, the bubble adapts its shape to the IP electron probability density before 
decaying by photoemission to the deformed IS state. Consequently, the bubble configuration 
at the emission time can be obtained by minimizing the grandpotential of the system keeping 
the electron in the excited IP state. We have done it by solving Eqs. (ISHlj) taking for $ the 
$iP wavefunction. In this case, a Gram-Schmidt scheme has been implemented to determine 
both the IP and IS relaxed states that obviously no longer correspond to a spherical bubble. 
In this axially symmetric environment, the spherical nL states are split depending on the 
value (m) of the orbital angular momentum on the symmetry z-axis, and the ±m states are 
degenerate. We have found that, within a nL manifold, the m states are ordered in increasing 
\m\ values. For this reason, we will refer to the axially symmetric state that corresponds 
to the m = submanifold when we speak of a deformed 'nL' state. When needed, we shall 
use the conventional notation for the orbital angular momentum of single particle states 
in linear molecules, namely a, tt, 0, . . . for \m\ = 0, 1, 2, 3, . . ., and superscripts + (— ) for 
specularly symmetric (antisymmetric) states. 

Figure [H displays quasi-equilibrium ebubble configurations at different stages of the 
absorption-emission cycle obtained at P = 0. The electron probability densities are rep- 
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resented by colored clouds, these with one lobe correspond to IS states (spherical bubble, 
picture 1; deformed bubble, picture 4), and these with two lobes correspond to IP states 
(spherical bubble, picture 2; deformed bubble, picture 3). In this figure, the line indicates the 
bubble dividing surface, i.e., the surface at which the liquid density equals half the saturation 
density value po, e.g. 0.0218 at P = 0, and represents the surface of the helium bubble. 
We have found^^ that at P = 0, the IP^IS emission energy is 36 meV, close to the 35 
meV found in Ref. id, constituting another excellent test of the theoretical framework used 
by us and by these authors. The energy released in the optical IS- IP absorption-emission 
cycle can be determined combining the results we have obtained in previous works.-^^^ For 
instance, at P = the released energy is 69 meV (compare with the 76 meV obtained in 



Ref. 



id). This energy is transferred to the superfiuid through generation of different kind 



of excitations. 

Quasi-equilibrium configurations of the ebubble relaxed around the IP state are shown 
in Fig. [2] for several P values. In this figure, helium is represented by warm colors, and the 
electron probability density (arbitrary units) by cool colors. The relaxation of the bubble 
around the IP state produces a characteristic two-lobe peanut structure whose waist -or 
neck- is progressively marked as the pressure applied to the liquid increases. Notice that 
helium displays a stratified density around the bubble. This feature appears whenever the 
superfiuid presents a kind of free surface, as in drops, films or bubbles 

Figure [2] shows that at a pressure of ^ 8 bar helium starts to penetrate between the two 
lobes of the electron wavefunction. At ^ 9 bar the helium density in this region reaches 
the saturation density, and the bubble splits into two baby bubbles. This produces an 
abrupt change in the emission energy, falling an order of magnitute between P = 8 bar 
and P = 9 bar; in the 'broken neck' region extending up to the solidification pressure, the 
photon emission energy is barely ^ 1 meV.^^ This is expectable at these pressures, as the 
main difference between the IS and IP probability densities appears in the waist region. If 
this region is inaccessible to the electron due to the presence of helium, these states become 
almost degenerate. On the contrary, if this region is not accesible to the superfiuid due, 
e.g., to the presence of a vortex whose vorticity line coincides with the symmetry axis of the 
ebubble, the baby bubbles may be held together by a tiny neck.^^ 

It is worth pointing out that some of the quasi-equilibrium configurations displayed in 
Fig. [2] may not be reachable in the evolution of the bubble. The reason is that helium falling 
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in the waist region during the violent coUapse may produce a large pileup of superfluid in 
that region, thus causing the actual breaking of the neck at pressures much smaller than the 
8 bar obtained quasi-statically. This possibility has been anticipated by Maris .- 



III. TIME EVOLUTION OF THE EBUBBLE: PICOSECOND DYNAMICS 

The dynamics of the excess electron localization in liquid helium has been adiabatically 
addressed by Rosenblit and Jortner using a sharp surface model for the bubble. ^''^^^ The 
superfluid was considered as incompressible, and the bubble expansion time, i.e., the time 
for creating the ebubble, was estimated to be 8.5 ps when energy dissipation by emission of 
sound waves was taken into account. This process exhibits a marked P dependence, the 
higher the pressure the shorter the expansion time. 

Later on, the incompressibility approximation was relaxed using a DF approach, flnding 
that at P = the bubble surface breathes with a period of about 130 ps.^ These calcu- 
lations have revealed that the localization process may launch shock waves, and that the 
subsequent main dissipation mechanism is sound radiation; excitations in the roton well 
were not produced.^ In the present work we consider that the ebubble has had time enough 
to relax to its spherical IS ground state and the electron is subsequently excited by light 
absorption to the IP state, whose dynamical evolution is the subject matter of this section. 



A. Adiabatic time evolution 



Since the electron evolves much faster than helium as their mass ratio is rriHe/'^e — 7300, 
we have followed the dynamics of the excited ebubble by combining the actual time evolution 
of the liquid with an adiabatic evolution for the electron. Within this approximation, the 
electron wavefunction is found by solving, at each time step, the static Schrodinger Eq. 
(II]) to obtain the instantaneous IP electron state $ip, and the structure of the liquid is 
obtained by determining the complex, time-dependent effective wavefunction ^^(r,t) from 
the time-dependent DF equation 

ot 2 rriHe 

+ |w[p,v] + |$ipr^^:^^U(r,t) , (5) 
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where the effective potential v] is given e.g., in Refs. l30l]3ll and has an exphcit depen- 
dence on the local current field j(r) = p(r) v(r) arising from the backfiow term the static 
potential energy in Eq. ([T]) lacks of. These coupled equations are solved imposing as initial 
conditions the stationary solution of the superfiuid for the spherical IS ebubble, employing 
the electron wavefunction $ip for building the e-He interacion, as indicated in Eq. ([5j). The 
time step has been set to 0.01 ps, and we have used a fourth order Runge-Kutta method 
to obtain the first time steps of the evolution. To determine the solution for subsequent 
times, we have used Hamming's (predictor-modifier-corrector) method. This procedure is 
very robust and accurate, even for large amplitude motions.— 

B. Results 

We have solved the adiabatic-dynamic coupled equations for P = 0, 0.5, 1, 2, 3, and 5 
bars. We will mostly show results for the two extreme pressure values, namely and 5 bars. 
The evolution starts by stretching the bubble along the symmetry z-axis, and shrinking its 
waist. This produces density waves in the liquid that take away a sizeable part of the energy 
injected into the system during the absorption process, 105 meV at P = and 148 meV at 
5 bar.- 

The evolution can be safely followed for about 25-30 ps. For larger times, the density 
waves refiected on the box boundaries get back to the region where the bubble sits, spoiling 
the calculation. This time interval is large enough to see bubble splitting at the higher pres- 
sures. Otherwise, one needs to introduce a source of damping in the equation governing the 
liquid evolution [Eq. ([5])], to prevent sound waves from bouncing back. Usually, introducing 
damping requires to enlarge the calculation box to accomodate a buffer region where waves 
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This increases the number of grid points and 



are washed out, see e.g. Refs. Is,! 
slows the calculation. 

Below 1 bar, we have found that the ebubble configuration is simply connected and ra- 
diates a sizeable part of the excitation energy as sound waves. For instance, at P = 
bar, the energy difference between the spherical IP configuration and the relaxed IP quasi- 
equilibrium configuration is ^ 40 meV (see Fig. [1]). The ebubble undergoes damped oscil- 
lations that will likely lead it to the corresponding quasi-static IP configuration described 
in Sec. II. As a consequence, it would eventually decay radiatively to the deformed IS state 
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that will radiationless evolve towards the spherical IS state. 

An example of this sort of evolution is shown in Fig. [3] for P = 0. We have found that 
after 15 ps, the shape of the IP bubble is similar to the quasi-static configuration referred 
to in Sec. II. Using a simpler model, Maris has found a smaller value, 11 ps.^ The difference 
is a natural consequence of the two basic approximations he has made, namely treating the 
liquid as incompressible and neglecting sound wave radiation. Due to the inertia of the 
bubble in the expansion process, it continues to stretch in the direction of the symmetry 
axis. This dilatation in the z-direction goes on for the largest times we have followed the 
evolution (30 ps), accompanied by the appearance of a more marked neck. 

At 1 bar, the neck collapses due to the large kinetic energy of the liquid filling in the 
region between the two IP lobes, and the ebubble configuration becomes non-simply con- 
nected. This causes the -deformed- IP and IS levels to become nearly degenerate, and their 
probability densities are almost identical. The appearance of any asymmetric fiuctuation, 
which is beyond the scope and capabilities of our framework, will cause the electron to even- 
tually localize in either of the baby bubbles. The subsequent evolution of the system is the 
collapse of the empty baby bubble and the evolution of the other one towards the spherical 
IS ground state. In this case, the excited IP bubble decays to the IS spherical configuration 
without passing through the IP quasi-static configuration described in the previous section, 
and the de-excitation is non-radiative. An example of this sort of evolution is shown in 
Fig. [4] for P = 5 bar. For this pressure, we have found that a configuration similar to the 
simply connected quasi-static one is attained after 10 ps, although the inertia of the bubble 
expansion breaks the quasi-static neck at about 18 ps. The density pileup in the neck region 
continues and at about 22 ps the helium density in this region has a peak of ^ 2po, whose 
relaxation pushes the two baby bubbles in opposite directions helping the fission process. 

Our calculations are in agreement with cavitation experiments^^ indicating that relaxed 
IP bubbles are only produced for pressures smaller than about 1 bar, and that above this 
value the decay of the excited bubble has to proceed differently, likely radiationless. Indeed, 
we have found that the IP bubble fissions at P = 1 bar, but it does not at P = 0.5 
bar. Our results are also in agreement with the interpretation^'' of the vanishing of the 
photoconductivity signal below 1 atm experimentally observed by Grimes and Adams.- 
According to this interpretation, an ebubble in the IP state is unstable against a radiationless 
de-excitation back to the ground state, the electron ultimately settling into one of the baby 
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bubbles while the other collapses with phonon ('heat') emission. It is this relased heat that 
drives the photocurrent. Below that pressure, the ebubble decays radiatively, it does not 
release enough heat, and is not detected in the photocurrent experiment. 

The evolution of the electron energies for the IS (empty) and IP (occupied) states, 
together with a representation of the electron probability densities, is presented in Fig. [5] 
for two pressure values. 

The fission of the bubble at P = 1 bar happens after 60 ps. To obtain this result, we 



have proceeded as in Refs. 
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221 . introducing a damping term in Eq. ([5j). To make sure 



that the bubble does not fission at P = 0.5 bar, we have also introduced a damping term 
for this pressure. 

We have studied the excitations produced in the liquid by the expansion of the ebubble. 
From the evolution of the first wave front, we have estimated that it moves in the z direction 
at ^ 330 m/s at P = 0, and at ~ 410 m/s at P = 5 bar. These values are well above the 
speed of sound in helium at these pressures, meaning that the dynamics is highly non-linear. 
Besides, we have determined the nature of these excitations by Fourier analyzing the density 
waves along the z axis, leaving aside the region near the bubble surface. The density profile 
is shown in the top panel of Fig. [6l corresponding to P = bar at 13 ps, and the Fourier 
transform of the density fiuctuation (related to the static structure function of the liquid) 
is presented in the bottom panel. Apart from the low-g component, arising from the mean 
density profile, one can identify two distinct peaks, the more marked one at g ^ 0.8 A-i m 
the phonon branch near to the maxon region, and another at g ^ 2.3 slightly to the 
right of the roton minimum. A similar 'roton' peak was found in Ref. Is'. A less marked 
peak appears at g ^ 1.7 A~^, slightly to the left of the roton minimum. Similar peaks have 
been found for shorter and larger times. From the relative intensity of these peaks, we are 
prone to identify most of the emitted waves as high energy 'phonons'. 

We have also analyzed the effect of the backfiow term on the appearance of the density 
waves. One can see from the bottom panel of Fig. [6]that neglecting this term changes a little 
the relative intensity of the phonon and roton peaks, increasing the former and decreasing 
the later as expected from the eflFect of the backfiow term on the excitation modes of the 

We want to stress that rotons are not excited if one uses 

In this case, only the low-g 



superfiuid, see Fig. 12 of Ref. 



17. 



a less accurate, local functional such as that of Refs. 



18 
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phonon spectrum of the superfiuid is quantitatively reproduced. Whether this has a sizeable 
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influence on the ebubble dynamics or not, can only be ascertained by a detailed comparison 
between the results obtained using both kind of functionals. 

IP ebubbles may be excited by photoabsorption to the ID state, either to its m = 
component (2(7+) or to its m = ±1 components (Ivr"). The absorption spectrum at different 
stages of the time evolution can be measured in a pump-probe experiment by which the 
ebubble is excited by two consecutive laser pulses. The delay set between these pulses should 
correspond to the time elapsed between the excitation of the spherical bubble and the time 
at which the absorption spectrum of the IP ebubble is recorded. Time-resolved electronic 
spectroscopy has been proposed by Rosenblit and Jortner as a tool for the exploration of 
the localization dynamics of the excess electron.— 

Time-resolved excitation energies are shown in Fig. [3 at P = and 5 bar. While 
the evolution of the la~ Itt" excitation is qualitatively similar at both pressures, the 
la~ 2(7+ excitation evolves differently in the high pressure regime when the bubble 
splits. Indeed, at zero bar the excitation energy smoothly decreases with time, whereas at 5 
bar it starts decreasing, increasing next, and eventually becoming larger than the excitation 
energy to the Itt" state. Note that both the change in behaviour and the crossing take place 
some picoseconds before bubble splitting. 

The 1(7~ Itt" transition is little affected by bubble splitting because it involves two 
states with negative specular symmetry, which means that the probability density of both 
states is zero in the neck region. On the contrary, the la~ 2(7+ transition involves two 
states with different specular symmetry and thus it is more affected by bubble splitting. 

Along with the excitation energies, some values of the associated oscillator strengths are 
displayed in Fig. [71 The oscillator strengths have been calculated in the dipole approxima- 
tion as: — 

fa, = ^{Ea-E,)\{a\r\b)\\ 

As known, the oscillator strengths fulfill a sum rule that in the one-electron case is fab = 
1.— At both pressures, we have found that these transitions have comparable oscillator 
strengths. The largest difference appears for P = 5 bar in the split-bubble regime. In it, the 
strength of the la~ 2(7+ transition is roughly half that of the la~ Itt" transition. We 
thus conclude that the analysis of the peak energy and oscillator strength of the la~ 2(7+ 
transition might disclose the fission-like behaviour of the excited IP bubble, complementing 
the experimental information gathered from cavitation and photoconductivity experiments. 

11 



The current field [j(r) = p(r)v(r)] is sliown in Fig. [8] for P = and 5 bar at 12 and 
22 picoseconds. At 12 ps the current fields are qualitatively similar for both pressures: the 
bubble expands along the symmetry axis and shrinks in a plane perpendicular to it. At 22 
ps, when P = 5 bar, large currents keep bringing liquid into the neck region, splitting the 
bubble and producing important density oscillations in the central region. 

We have also followed the collapse of the 2P bubble at P = 0. About 2 ps after the 
collapse has started, a sizeable part of the excitation energy has been released into the 
liquid and two waves are distinguishable around the bubble, as shown in Fig. M These 
waves travel through the liquid at the same speed as in the IP bubble case, ^ 330 m/s. 
Shortly after 7 ps, the m = levels of the 2P and IF states become very close and, as 
discussed in the next section, see Eq. ([7j), the adiabatic approximation fails. At this point, 
the 2P bubble displays an incipient four-lobe shape arising from a similar structure in the 
2P electron probability density. It is worth mentioning that a likely related effect, namely 
the n^ar degeneracy of the 2P and IF states, was found in the quasi-static calculations of 
Ref. I26I as P increased. 

In view of the mentioned failure and the lacking of experimental information on the 
collapse of the 2P bubble, we have closed its study at this point, leaving it for future work. 

IV. VALIDITY OF THE ADIABATIC APPROXIMATION 

The validity of the adiabatic approximation in the first stages of the bubble collapse, 
when the topology of the bubble is simply connected, stems from the very different time 
scale of the electron motion as compared to that of the bubble, represented by the period 
of its shape oscillations. If the fluid is incompressible and the bubble spherical, the surface 
A-modes of the cavity are at energies 



where 7 and po are the surface tension and atom density of the liquid, respectively. For 
A = 2 this energy is about 1 K, and the period of the oscillation is ^ 50 ps. 

The situation may change in the course of the collapse because the energy difference 
A = Eip — Eis between the deformed states decreases and the time scale r = h/A may 
become similar to the period of the shape oscillations of the deformed bubble. Since A is 
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small in the two bubble regime, the approximation likely fails there.— »^ It is worthwhile 
mentioning that neck fluctuations, not included in our approach nor in previous works, 
would pinch off the bubble at earlier stages of the collapse, in a similar way as they may 
cause the prompt scission of the flssioning atomic nucleus after the saddle conflguration has 
been overcome.— 

On the light of our model, in which no assumptions are made on the shape of the bubble 
nor the impenetrability of the bubble surface by the localized electron, it is instructive to 
analyze the validity of the adiabatic approximation assuming that, during the evolution, the 
bubble keeps its original axial symmetry and specular symmetry about the plane perpendic- 
ular to the symmetry axis that contains the node of the IP state. This excludes any possible 
fluctuation and the appearance of asymmetric modes, like the breathing mode discussed in 
Ref. Our discussion relies on the detailed presentation by Messiah, that we summarize 
in the following. 

The subsystem to which the adiabatic approximation is applied is the electron, whose 
wavefunction is decoupled from that of the liquid. This wavefunction evolves in the potential 
fleld generated by the liquid distribution, and its Hamiltonian is time-dependent, 1-Le{t) = 
T-Le[pHe{t)]- Let 0^ be an eigenfunction of the Hamiltonian at time t, so that l-Le{t)4L = 
£n{t)(/)l^. If (pnit) is the actual wavefunction 0^ evolved up to time t, one has (/)n{t) = U{t)(/)^^ 
where U is the evolution operator. In the adiabatic approximation, one identifles (/>n(i) with 
0^, the intuitive justiflcation being that if one perturbes the subsystem slowly and gently 
enough, it has enough time to adapt itself to the new environment 'with no inertia' from 
the past conflguration. 

The error made in this approximation for a given state \i) -the validity of the adiabatic 
approximation is assessed for a given state of the subsystem, not necessarily for them all- is 
deflned as the probability of flnding the subsystem in a state different from the initial one 
evolved in time within the 'true' dynamics, rjij = This error can be written 

in a workable way as^- 

2 



I) 



(7) 



If Tjij <^l\/i ^ j, the adiabatic approximation is justifled. It is usually understood that it 
breaks down when the levels get very close, or when they cross each other. Notice however 
that this assumes that these states can be connected by the evolved Hamiltonian. If a 

13 



symmetry is dynamically conserved (in our case, angular momentum and specular symmetry 
are), then the adiabatic evolution of states with a given quantum number associated to this 
symmetry is not perturbed by states with different values of this quantum number. Although 
sometimes ignored, this is a very reasonable statement. 

In the case of the la~ state arising from the spherical IP manifold, the adiabatic ap- 
proximation holds even when its energy becomes almost identical to the energy of the la^ 
state arising from the spherical IS one, i.e. a small A does not invalidate the adiabatic 
approximation. The closest state having the same angular momentum and specular sym- 
metry is the 2a~ one arising from the spherical IF manifold. At P = 5 bar, we have found 
that these two states are 2000 K apart in the 10-20 ps range. Since h ^ 7.6 Kps, one has 
h/{E2cr- — Ei^-Y ^ 2 X 10~^ psK~-^. The value of the matrix element in Eq. ([7j) is some 
tens of kelvin per picosecond, so that the adiabatic approximation would be fulfilled even 
for the configuration displayed in Fig. [Hat 25 ps. Indeed, we have calculated ij in the above 
time range and have found that it is of the order of 10~^. 

It is also worth analyzing the stability of the quasi-equilibrium configurations when the 
symmetries are not exactly conserved because of perturbations from the environment. In 
this situation, let us assume that when the bubble splits the electron localizes in one of the 
lobes. Leaving out the discussion on the actual localization process, we have tried to infer 
the likely evolution of an ebubble with a localized electron. The localized electron state in 
either baby bubble is approximated by: 



Consider now a short-time dynamics in which the liquid is kept frozen. The evolution of, 
e.g., the localized state is an oscillation between the two lobes 



where u) = {Eip — Eis)/h. If this frequency is large enough, the liquid cannot react to the 
localization of the electron in either lobe, and will essentially behave as if the electron were 
delocalized. 

The time elapsed between two consecutive localizations of the electron in a given baby 
bubble isT = tt/uj. The value of this period as a function of pressure for the quasi-equilibrium 
configurations is displayed in Fig. [101 In the split-bubble regime (P ^ 9 bar), this period is 




(8) 



(9) 
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of several ps, indicating that the electron localization dynamics into one of the baby bubbles 
is not a trivial process to address. The electron will bounce back and forth as the liquid 
tries to adapt to it. Real time calculations are thus needed to describe electron localization. 

It is clear that the previous discussion on the validity of the adiabatic evolution lacks for 
incorporating fluctuations or excitations of low energy modes that may appear in the course 
of the bubble evolution and couple the IP and IS states that otherwise are not, as previously 
discussed. One such mode has been thoroughly addressed by Elser:^'' a peanut conflguration, 
whose walls are impenetrable by the electron, is represented by two intersecting sharp spheres 
of radius R2 (instead of the deformed baby bubbles displayed in Fig. Hj) joined along a 
circular oriflce of radius a. These spheres are breathing in counterphase, producing an 
antisymmetric mode whose stiffness JC and inertia Ai can be obtained analytically. This 
mode is very appealing, as it represents a small, swifting imbalance of the symmetric electron 
probability density. 

In the harmonic limit, if a <C i?2 5 the stiflFness and inertia of the asymmetric mode are— 



K = 16 



2 



:)3 



M 47r^m4HePo^2 • (10) 

In these equations, Eq = jR^ ^ 9 meV represents the energy unit, with 7 = 2.36 x 10~^ 
meVA~^ being the surface tension of the liquid, i?o ^ 20 A is the radius of the spherical 
bubble, i?2 ^ 16 A is the radius of the baby bubbles, ^ 1.70 is a dimensionless constant, 
and A = Eip — Eis. The frequency of the antisymmetric breathing mode is given by 



^AB = y/C/A^, and the radius of the oriflce is^'' 

The adiabatic approximation fails when ujab — A//i. This yields A ^ 0.14 meV in the 
P = — 5 bar range, as only the flrst term in the stiffness turns out to be relevant in this 
regime. Thus, keeping only the flrst term in /C, one gets^^ 

1/9 

-0.13 . (12) 



a 

'R2 



( 27 \ ( me \ /Ra\ 



Hence, a ~ 2.1 A. Clearly, such analytical results cannot be obtained within the DF ap- 
proach, but we can use them to determine whether the dynamic and static configurations 
shown in Sees. II and III are reliable. 
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The adiabatic approximation thus holds at P = 0, as the neck radius is fairly large, see 
Fig. [31 and A is always much larger than 0.14 meV, see Fig. [51 From Fig. [H we also 
conclude that, at P = 5 bar, the adiabatic approximation is valid up to nearly the collapse 
of the waist. Indeed, the neck radius of the helium configuration at about 17-18 ps is ^ 2 
A, see Fig. [H It is worth noting the difficulty in defining an effective radius for the orifice 
when the surface of the bubble is diffuse; we recall that the surface thickness of a ^He drop 
of 10^-10^ atoms is some 6-8 A.^^ Note also that the surface thickness of the helium bubble 
is rather independent of the curvature of the surface, as can be inferred from the fairly 
constant bright region around the bubbles displayed in Figs. [3]and[ll 

Since we do not treat the bubble as impenetrable to the excess electron, the relation 
between the actual A and a values should not exactly be as given by Eq. ( fTTI) . Using 
the result hujAB < A = 0.14 meV as a criterion for the applicability of the adiabatic 
approximation instead of reaching the limiting value a ^ 2 A, we find that the approximation 
holds up to 21 ps, when the bubble has already split into two baby bubbles. Both procedures 
indicate that when the adiabatic approximation likely fails, the baby bubbles have already 
developped. 

The previous analysis leads us to conclude that, at high pressures, baby bubbles are 
formed some tens of picoseconds after the starting of the collapse of the IP bubble. From 
this point on, the likely fate of the system is the localization of the electron in one of the baby 
bubbles and the collapse of the other. This process is helped/triggered by fiuctuations that 
break the specular symmetry of the ebubble configuration. As mentioned, determining the 
time scale of electron localization is beyond the capabilities of the adiabatic approximation. 
It has been calculated^^ that once the electron is localized, it takes to the superfiuid some 
20 ps to adapt to it while the other baby bubble is absorbed. 

V. SUMMARY 

Within density functional theory, we have carried out an analysis of the adiabatic evolu- 
tion of the excited electron bubble in superfiuid liquid ^He. We have found that for pressures 
below 1 bar, the IP ebubble may relax to its quasi-static equilibrium configuration and even- 
tually decay radiatively to the deformed IS state. This state evolves non-radiatively to the 
spherical IS bubble, completing the absorption/emission cycle. This conclusion arises in 
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part from studies carried out for one hundred picoseconds using a less accurate functional,— 
whose results qualitatively agree with ours for the first tens of picoseconds. 

At higher pressures, the situation drastically changes and the excited IP bubble no longer 
decays to the quasi-static equilibrium configuration, whose physical realization is unlikely. 
Indeed, our analysis of the adiabatic approximation indicates that it is valid up to a point 
where two deformed, nearly disconnected baby bubbles appear in the dynamical evolution, 
pointing towards a fission-like de-excitation process, the likely subsequent evolution of the 
system being the localization of the electron in one of the baby bubbles and the collapse of the 
other. This collapse takes some 20 ps,^^ and the whole de-excitation process is radiationless. 

We have also found a marked change in the behavior of the time-resolved absorption spec- 
trum of the IP bubble depending on whether the bubble fissions or not, i.e., on the liquid 
pressure. This change is in principle an experimentally accessible observable whose deter- 
mination may complement the information obtained from cavitation and photoconductivity 
experiments. 

Our analysis of the collapse of the 2P bubble has shown that the adiabatic approximation 
breaks down at an early stage of the dynamical process due to the crossing of the 2P and 
IF states. Although disclosed by the adiabatic approximation, this crossing has nothing to 
do with the approximation itself, but is inherent to the dynamics of the electron bubble. 
From the crossing point on, the bubble will relax around a mixed state with 2P and IF 
components, and hence the physical realization of a pure quasi-equilibrium 2P configuration 
is unlikely. It is very plausible that the same applies to other high energy 'nL' ebubbles 
generated in the absorption process. The possibility that some of them undergo a sponta- 
neous symmetry breaking, as suggested by Grinfeld and Kojima for the 2S state, can only 
reinforce our conclusion. Obviously, this does not question the existence of either relaxed 
quasi-equilibrium configurations at low pressures, or of baby bubbles at high pressures, aris- 
ing from the evolution of the spherical 2P bubble. It just means that, on the one hand, the 
relaxed bubble will not be a pure 2P configuration and, on the other hand, to study the 
de-excitation of these bubbles one has to go beyond the adiabatic approximation and carry 
out a more demanding real time dynamics calculation for the electron. 



17 



Acknowledgments 



We thank Humphrey Maris, Francesco Ancilotto, Dafei Jin, and Alberto Hernando for 
useful discussions. This work has been performed under Grants No. FIS2008-00421/FIS 
from DGI, Spain (FEDER), and 2009SGR1289 from Generalitat de Catalunya. 



1 M. Rosenblit and J. Jortner, J. Chem. Phys. 124, 194505 (2006); ibid. 124, 194506 (2006). 

2 H. J. Maris, J. Phys. Soc. of Jpn. 77, 1 (2008). 

^ J. A. Northby and T. M. Sanders, Phys. Rev. Lett. 18, 1184 (1967). 
^ W. B. Fowler and D. L. Dexter, Phys. Rev. 176, 337 (1968). 

5 T. Miyakawa and D. L. Dexter, Phys. Rev. A 1, 513 (1970). 

6 C. C. Grimes and G. Adams, Phys. Rev. B 41, 6366 (1990). 
^ C. C. Grimes and G. Adams, Phys. Rev. B 45, 2305 (1992). 

s J. Eloranta and V. A. Apkarian, J. of Chem. Phys. 117, 10139 (2002). 

9 V. Grau, M. Barranco, R. Mayol, and M. Pi, Phys. Rev. B 73, 064502 (2006). 

L. Lehtovaara and J. Eloranta, J. Low. Temp. Phys. 148, 43 (2007). 

^ B. DuVah and V. Celh, Phys. Rev. 180, 276 (1969). 

2 M. Pi, R. Mayol, A. Hernando, M. Barranco, and F. Ancilotto, J. Chem. Phys. 126, 244502 
(2007). 

^ W. Guo, D. Jin, G. M. Seidel, and H. J. Maris, Phys. Rev. B 79, 054515 (2009). 

A. Ghosh and H. J. Maris, Phys. Rev. Lett. 95, 265301 (2005). 
^ J. Tempere, L F. Silvera, and J. T. Devreese, Sur. Sci. Rep. 62, 159 (2007) 
^ J. Fang, A. E. Dementyev, J. Tempere, and L F. Silvera, Rev. Sci. Instrum. 80, 043901 (2009) 
^ F. Dalfovo, A. Lastri, L. Pricaupenko, S. Stringari, and J. Treiner, Phys. Rev. B 52, 1193 

(1995). 

^ S. Stringari and J. Treiner, Phys. Rev. B 36, 16 (1987); J. Chem. Phys. 87, 5021 (1987). 
9 N. G. Berloff and P. H. Roberts, Phys. Rev. B 63, 024510 (2000); J. Phys. A 34, 81 (2001) 
'0 E. J. M. Madarassy and C. F. Barenghi, J. Low Temp. Phys. 152, 122 (2008). 

D. Jin, W. Guo, W. Wei, and H. J. Maris, J. Low Temp. Phys. 158, 307 (2010). 

D. Jin and H. J. Maris, J. Low Temp. Phys. 158, 317 (2010). 



18 



23 N. G. Berloff and P. H. Roberts, J. Phys. A 32, 5611 (1999). 

24 N. G. Berloff and P. H. Roberts, Phys. Lett. A 274, 69 (2000) 

25 D. Mateo, A. Hernando, M. Barranco, and M. Pi, J. Low Temp. Phys. 158, 397 (2010). 
2*^ H. J. Maris, J. Low Temp. Phys. 132, 77 (2003). 

27 V. Elser, J. Low Temp. Phys. 123, 7 (2001). 

28 R. Jackiw, C. Rebbi, and J. R. Schrieffer, J. Low Temp. Phys. 122, 587 (2001). 

29 F. Ancilotto, M. Barranco, and M. Pi, Phys. Rev. Lett. 91, 105302 (2003). 

30 L. Lehtovaara, T. Kiljunen, and J. Eloranta, J. of Comp. Phys. 194, 78 (2004). 

31 L. Giacomazzi, F. Toigo, and F. Ancilotto, Phys. Rev. B 67, 104501 (2003). 

32 F. Ancilotto, F. Dalfovo, L. P. Pitaevskii, and F. Toigo, Phys. Rev. B 71, 104530 (2005). 

33 E. Cheng, M. W. Cole, and M. H. Cohen, Phys. Rev. B 50, 1136 (1994); Erratum ibid. 50, 
16134 (1994). 

34 M. Pi, M. Barranco, R. Mayol, and V. Grau, J. Low Temp. Phys. 139, 397 (2005). 

35 M. Prigo and S. G. Johnson, Proc. IEEE 93, 216 (2005). 

36 M. Barranco, M. Guilleumas, E. S. Hernandez, R. Mayol, M. Pi, and L. Szybisz, Phys. Rev. B 
68, 024515 (2003). 

37 M. Rosenblit and J. Jortner, Phys. Rev. Lett. 75, 4079 (1995). 

38 M. Rosenblit and J. Jortner, J. Phys. Chem. A 101, 751 (1997). 

3^ A. Ralston and H. S. Will, Mathematical methods for digital computers (John Wiley and Sons, 
New York, 1960). 

M. Barranco, R. Guardiola, S. Hernandez, R. Mayol, J. Navarro, and M. Pi. J. Low Temp. 
Phys. 142, 1 (2006). 

41 Ch. Cerjan, D. Kosloff, R. Kosloff, and M. Reshef, Geophysics 50, 705 (1985). 

42 H. J. Maris, A. Ghosh, D. Konstantinov, and M. Hirsch, J. Low Temp. Phys. 134, 227 (2004). 

43 M. Weissbluth, Atoms and Molecules (Academic Press, New York, 1978). 

44 Recall that if the initial state |6) is not the ground state, the sum includes negative terms arising 
from de-excitation transitions, hence a partial sum can be greater than unity. 

45 H. J. Maris, J. Low Temp. Phys. 120, 173 (2000). 

46 A. L M. Rae and W. F. Vinen, J. Low Temp. Phys. 123, 1 (2001). 

47 U. Brosa, S. Grossmann, A. Miiller, and E. Becker, Nucl. Phys. A 502, 423 (1989). 

48 A. Messiah, Quantum Mechanics (North-Holland, Amsterdam 1962). 

19 



271 does not conserve the volume of the bubble, the stiffness 



271 if P 7^ 0. In the numerical application we 



Since the asymmetric mode of Ref. 
slightly changes from the expression given in Ref. 
have used the 7(P) law deduced in Ref. 

This equation differs from Eq. (9) of Ref. Q because of an error made there to obtain it from 
a straightforward combination of its Eqs. (3) and (8), as can be easily checked. 
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FIG. 1: (Color online) Ebubble quasi-equilibrium configurations at different stages of the 
absorption-emission cycle corresponding to P = 0. The electron probability densities are rep- 
resented by colored clouds. The dashed line is the bubble dividing surface. The size of the frames 
is 70 A x70 A. 
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FIG. 2: (Color online) IP bubble quasi-equilibrium configurations for P = 0, 5, 8, and 9 bar. 
Helium is represented by warm colors, and the electron probability density (arbitrary scale) by 
cool colors. The size of the samples is 70 A x70 A. 
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FIG. 3: (Color online) Adiabatic evolution of the IP ebubble at P = 0. The panels display the 
helium configurations corresponding to 0, 5, 10, 15, 20, and 25 ps. The dashed line represents the 
dividing surface of the quasi-equilibrium configuration at P = shown in Fig. [2l 
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FIG. 4: (Color online) Adiabatic evolution of the IP ebubble at P = 5 bar. The panels display 
the helium configurations corresponding to 0, 5, 10, 15, 20, and 25 ps. The dashed line represents 
the dividing surface of the quasi-equilibrium configuration at P = 5 bar shown in Fig. [2l 
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FIG. 5: (Color online) Adiabatic evolution of the energies of the IS (empty) and IP (occupied) 
states, together with a representation of the electron probability densities for P = (top panel) 
and P = 5 bar (bottom panel). The electron probability densities are also displayed at four selected 
time values. 
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FIG. 6: (Color online) Top panel: Density profile of the superfluid corresponding to P = bar at 
13 ps (only the z > part is shown). The inset shows the ebubble configuration in a 140 A X 140 
A frame, and the region displayed in the rectangle is the Fourier tranformed one. Bottom panel: 
Solid line, Fourier transform of the superfuid density profile shown in grey in the top panel. The 
dashed line is the Fourier transform of the density obtained at the same evolution time without 
including backflow effects. In both cases the peak at g = A~^ arises from the mean density 
profile, not subtracted from the local density before transforming. The lines have been drawn as a 
guide to the eye. 
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FIG. 7: (Color online) Time-resolved excitation energies at P = and 5 bar for electron transitions 
from the la~ state arising from the spherical IP, m = one, to either the Itt" state (circles), or 
to the 2cr^ one (squares), both initially belonging to the spherical ID manifold. In each panel, the 
oscillator strength for two selected time values is also displayed. Lines have been drawn as a guide 
to the eye. 
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FIG. 8: (Color online) Current field at P = (top panels) and 5 bar (bottom panels) at time 
t = 12 ps (left) and 22 ps (right). In each panel, the larger the arrow, the larger the current. 
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FIG. 9: (Color online) Adiabatic evolution of the 2P ebubble at P = 0. The panels display the 
helium configurations corresponding to 2, 3, 5, and 7 ps. 
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FIG. 10: (Color online) Top: superposition of the IS and IP states corresponding to the quasi- 
equilibrium configuration at P = 9 bar to approximately localize the electron in one of the 
baby bubbles. Bottom: tunneling period of a localized electron in one of the lobes of the quasi- 
equilibrium ebubble as a function of pressure. 
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